Dark Matter and Higgs Boson in a Model with Discrete Gauge Symmetry 
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In view of ongoing measurements of the Higgs-like boson at the LHC and direct searches for dark 
matter, we explore the possibility of accommodating the potential results in a simple new-physics 
model with discrete gauge symmetry as well as light neutrino masses. Specifically, we study collider 
and relic-density constraints on the new gauge coupling, predict the cross section of the dark matter 
scattering off nucleons, and compare it with current direct search data. We also discuss some of 
the implications if the dark matter is light. The new gauge sector of the model allows it to be 
compatible with the latest LHC information on the Higgs-like particle and simultaneously satisfy 
the requirements in its dark matter sector. 
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I. INTRODUCTION 

The particle with mass around 125 GcV recently dis- 
covered at the LHC has characteristics which, based on 
the measurements reported so far [l[ , suggest that it may 
be the Higgs boson of the standard model (SM). Whether 
that is indeed the case will likely become clearer in the 
near future after sufficiently more data are accumulated 
at the LHC. Needless to say, the imminent potential con- 
firmation on the Higgs' existence will have far-reaching 
implications for attempts to identify the nature of new 
physics beyond the SM. For it is widely accepted that new 
physics is necessary at least to explain the astronomical 
evidence for dark matter (DM) and the numerous exper- 
imental indications of neutrino mass 0. In anticipating 
the outcomes of the Higgs quest as well as upcoming re- 
sults of direct searches for DM, which are also ongoing, 
it is of great interest to explore a simple framework of 
new physics that can accommodate the possibilities. 

To do so, we adopt an economical model that offers not 
only a Higgs boson and DM of the popular weakly inter- 
acting massive particle (WIMP) type, but also a means 
to produce light neutrino masses. The parameter space 
of the model has ample room for a Higgs with SM char- 
acteristics or one that is less SM-like. The Higgs and DM 
sectors are linked in that the Higgs may decay substan- 
tially into DM particles if kinematically allowed. This 
makes the Higgs and DM searches complementary for 
probing the model. Another of its salient features is that 
the DM stability is realized under a discrete Z2 symme- 
try that is not imposed in an ad hoc manner, but instead 
emerges naturally as a remnant of an extra Abelian gauge 
group, U(l)^, spontaneously broken by the nonzero vac- 
uum expectation value (VEV) of a new scalar field 
We assume that the SM fermions carry U(l)^ charges, 
and this allows us to introduce right-handed neutrinos 
for gauge-anomaly cancellation and generating light neu- 



trino masses by means of the well-known seesaw mecha- 
nism [H, which is activated with the involvement of the 
same new scalar field. 

This paper is organized as follows. In Sec. [Hi we de- 
scribe our model in greater detail. We also briefly com- 
pare the important features of our model with a number 
of other scenarios in the literature. In Scc. lIIIi we ob- 
tain constraints on the new gauge sector from collider 
experiments. In Sec. lIVi we address further constraints 
on the model from the relic abundance data and predict 
the DM-nucleon scattering cross-section subject to direct 
detection measurements. We also discuss the invisible de- 
cay of the Higgs boson if the DM is sufficiently light and 
consider some of the implications in relation to the Higgs 
hunt at the LHC. Finally, Scc.lVl contains the summary 
of our work and some additional discussions. 



II. MODEL DESCRIPTION 



We introduce a minimal number of new particles for 



our purposes: three right-handed neutrinos, v kR , two 
complex scalar fields, S and D, which are singlets un- 
der the SM gauge group, and a Z' boson for the U(l)^ 
gauge symmetry. We collect the quantum number as- 
signments for the fcrmion and scalar fields of our model 
in TablelH One can therefore see that U(l)^ is none other 
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TABLE I: Charge assignments of the fermions and scalars in 
the model. / SM (<?g M ) denotes SM fermions (their assign- 
ments) and H the usual complex doublet. For quarks and 
leptons, Q = 1/3 and —1, respectively. 
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than U(1)b_l 0,0, with B and L referring to baryon 
and lcpton numbers, respectively. The kinetic terms of 5* 
and D along with a renormalizable potential V for them 
and H are given by 

C = (VDyV^D + ^SyV^S - V , (1) 
V = fi 2 D \D\ 2 - »l\S\ 2 + fx DS (D 2 Si +K.c) 

+ 2X DS \D\ 2 \S\ 2 + 2{X DH \D\ 2 + X HS \S\ 2 )H^H 
+ XolDf + X S \S\ 4 + (X H H^H - &)rfH , (2) 

where T>^ = + ig^Z'^ with the coupling and 
charge C associated with U(l)^. The H field develops 
a nonzero VEV as in the SM. The non-SM parts of V 
were discussed before {e.g., Refs. @,H). The role of DM 
is played by the D field — sometimes dubbed the darkon — 
via its lighter component. To maintain the Zi symmetry 
and hence the DM longevity, D must have zero VEV. The 
presence of both the S field and the /j, DS term is essential 
because the latter triggers the spontaneous breakdown 
U(l)<j — > Zi when S gets a nonzero VEV. 

The parameters in V should be chosen such that the 
vacuum has the above desired properties. Accordingly, 
we assume that all the A's in V are positive to render 
it bounded from below. Subsequently, upon expressing 
the VEV's of H and S as (H) = (0 v H ) T /V2 and 
(S) = v s /\/2, with v H s > 0, we arrive at 



V H(S) 



X 



X 2 



and so 
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^H(S) V H(S) + ^HS V S(H) > 



(3) 



(4) 



Furthermore, writing D = (DR+iDi)/y/2 in terms of its 
real and imaginary components leads to the combinations 



Dr>Dt 



(5) 



which are the squared masses of Dr and Dj. The mass 
difference between Dr and Z)j therefore depends on the 
sign of fi DS . We will take fi DS > so that Dj acts as 
the WIMP DM. The results of our analysis would be the 
same if we took /j,ds < 0, om y that the roles of Dr and 
Dj would be interchanged. 

After clcctroweak symmetry breaking, the remaining 
field h' in H = (0 v H + /i') T /%/2 will mix with s 1 = 

the mass eigenstates 



v s because of the X HS term in C. This results in 



h = h! cos 9 + s' sin 9 , ,s = s' cos 9 — h' sin 9 , (6) 
with the mixing angle 9 and masses m h s given by 



tan(2(9) = 



l HS 



M 2 H - MJ 



(7) 



2m? 8 = A^ + M| T [(Af| r -M|) 2 + M4 s ] 1/2 , (8) 



where 



M 



H,S 



= 2X 



H,S V H,S i 



Mfjs = 4A 



HS V H U S 



(9) 



The lighter state, h, is the physical Higgs boson. 

In the neutrino sector, the mass-generating terms have 
the form 



l \i v kR H Tt 2 L il - \^ki u 



KihR^mTS 1 + H.c. , (10) 



where k, I = 1,2,3 are summed over, r 2 is the second 
Pauli matrix, and Lil represents a lepton doublet. They 
give rise to the Dirac and Majorana mass (3x3) matrices 



M D = ^Xv H , M VR = ^X'v s , 



V2 ■ 



(11) 



respectively. Incidentally, the term in Eq. (|10|) can 
play the same role as the /i DS term in Eq. ^ for symme- 
try breaking. For the type-I seesaw mechanism to yield 
the light neutrino masses, the eigenvalues of Ai UR are ex- 
pected to be orders of magnitude bigger than a TeV Q • 
Thus, assuming that X H and X s in Eq. © arc roughly 
of similar order to the eigenvalues of A' in Eq. ([TP]) , we 
have M s > ITeV and also M s > M H . In addition, 
we will pick M HS -C M§ — M H , so that m h s ~ M H s 
and 9 -C 1, leading to h ~ h! and s ~ s'. It follows 
that we can neglect the effects of the heavy s on the pro- 
cesses of interest and the relevant interactions for Dr i 
are described by 

£d = -^ DH {D\ + D 2 )(v H h+\h 2 

~ \^d{D 2 r + D 2 ) 2 + g^D^Dj)^ . (12) 

From now on, we consider the possibility that Dr and 
Di are nearly degenerate. As a consequence, the DM 
relic abundance is determined not only by the annihila- 
tion rate of Dj, but also by that of Dr and/or both of 
them. In this so-called coannihilation case @ we assume 
specifically that 



(m Da 



m Dl )/r 



Dr 







(13) 



Therefore, 



Di(r)Dj( R ) 



the relevant reactions are mainly 
— > SM particles, from Higgs-exchange 
and contact diagrams, the latter if mDt.Dn > m h-> an< l 
DjDr Z'* — > SM fcrmions, as the fcrmions carry 
U(l)^ charges. Moreover, the X D part in C D is not 
pertinent to our purposes. 

Before proceeding to our numerical calculations, we 
would like to make a few remarks comparing this work 
to those in the literature. A number of earlier analyses 
addressed the phenomenology of some of the elements 
in our model separately, such as scalar DM in the ab- 
sence of a new gauge sector 0, Hol - [l3j or the Z' boson of 
a gauged U(1)b-l symmetry in conjunction with right- 
handed neutrinos, but no DM candidates Q. There are 
also studies dealing with models which possess DM sta- 
bilized by the remnant Z2 symmetry of a U(l) b-l that 
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is global [lH or local @, [TsJ . The Zi symmetry could in- 
stead be simplyput in by hand @, [iMl El , or it could 
be accidental [17|- Here we would like to emphasize that 
the model which we have adopted incorporates the min- 
imal mechanism for both stabilizing scalar DM with the 
Z-2. remnant of a gauged U(l) b-l and generating neu- 
trino mass through the spontaneous breaking of the same 
group. Furthermore, the most nontrivial aspects of our 
analysis — in distinction to others' — are that the Z 1 bo- 
son of the U(l) b-l as well as the Higgs boson contribute 
to the DM interactions with SM particles, and hence to 
DM annihilation, and that we treat the case where the 
Z' boson plays a dominant role in determining the DM 
relic density. As we will demonstrate later, this allows 
the model to provide not only a good candidate for DM 
that is consistent with the latest direct-search results, 
but also a Higgs boson that is SM-like in its couplings 
to the standard particles and, for sufficiently light DM, 
has a significant invisible decay mode compatible with 
the current LHC data. 



III. CONSTRAINTS ON THE NEW GAUGE 
COUPLING 



Since Xdh and in Eq. (fl"2"j) are free parameters, we 
will consider different interesting combinations of their 
contributions to the relic density and their potential im- 
plications for Higgs and DM direct searches. Before do- 
ing so, we first constrain g^ using other observables mea- 
sured at colliders. Since only S induces the Z' mass, 
m z , = 2g ( -v s , there is no tree- level Z-Z' mixing, im- 
plying that at tree level gQ has no effects on the Z-pole 
observables [TH, but can affect e + e~ and hadron colli- 
sions into fermion pairs. Neglecting Z-Z' kinetic mixing 
which can arise at loop level, we focus on e + e~ — > £ + £~ 
and the Drell-Yan (DY) process pp —> £ + £~X, which can 
restrict g^ well. 

Measurements of e + e~ — > £ + £~ for £ = /j,,t were per- 
formed at LEPII with center-of-mass energies from 130 
to 207 GeV [l9j]. We employ the data on the cross sec- 
tion and forward-backward asymmetry. Adopting their 
90% confidence-level (CL) ranges and employing the for- 
mulas give n in Ref. [18J, but with s-dependent Z and Z' 
widths [19j , we extract the upper limit on g^ as a function 
of m z ,, represented by the blue solid curve in Fig.[TJ 

The latest cross-section data on the DY process from 
the LHC [2(| reveal no deviation from the SM expecta- 
tions and hence no evidence of a Z' boson. We can derive 
an upper bound on the coupling constant g^ using the SM 
cross-section, following the method of Ref. [2l[ . The DY 
cross-section is numerically estimated using the CalcHEP 
package [HI by incorporating new Feynman rules in the 
model file. We count the events in the invariant-mass 
window of ±20% around the Z 1 mass for a luminosity 
of l.lfb -1 according to the recent experimental analy- 
sis [20| . This number of signal events is plugged into the 
one-bin log likelihood LL = 2[N \n(N / is) + is — N], where 
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FIG. 1: Upper limits on the gauge coupling of the Z 1 
boson versus its mass m z t from LEP II and LHC data on 



£ and pp 



v l X, respectively. 



N (is) is the number of events predicted by the SM (SM 
plus the Z' boson). A value of LL = 2.7 is taken for 
the 90% CL, and v is solved for. The upper limit on the 
cross section is then derived from the solved value of v 
for each Z' mass. This upper limit in turn constrains g^ 1 
as depicted by the purple dashed curve in Fig. Q] 



IV. DARK MATTER PHENOMENOLOGY 

In the simplest darkon model, with a real darkon and 
no new gauge sector, the m D < m h /2 region has mostly 
been disfavored by DM direct detection results as well 
as the recent observation of the Higgs-like resonance at 
the LHC, for the darkon-Higgs coupling is required to 
be sizable by the relic density data, rendering the Higgs 
mostly invisible [23|. In our model, the presence of the 
Z' boson alters the darkon phenomenology in impor- 
tant ways. Particularly, there are now Z'-mediated di- 
agrams, besides the Higgs-mediated ones, contributing 
to both darkon annihilation and darkon-nucleon interac- 
tions. One of the consequences is that the allowed pa- 
rameter space of the model can still comfortably make 
room for a light darkon, as we will sec below. 



A. Relic abundance 

In evaluating the darkon contributions to the relic den- 
sity Q D , we employ the relations |qL l2~H 



J = 



1.07 x 10 J 
y/(L~m pi J GeV 

~dx (<7eff ! rel) , 



In 



0.038 g efc m D m Pl (<T cS v rcl )(g*x f ) 1/2 , (14) 



where /i denotes the Hubble constant in units of 
lOOkm/s/Mpc, is the number of rclativistic degrees 
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of freedom below the freeze-out temperature Tj , m pl = 
1.22 x 10 19 GeV is the Planck mass, = m D /Tj with 
m D = rauj being the WIMP mass, (<J e fiV Ic i) = (av) is 
the thermally averaged product of the effective darkon- 
coannihilation cross-section and the relative speed of the 
darkon pair in their center-of-mass frame, and g eS is the 
darkon's effective number of degrees of freedom in the 
coannihilation case. Our choice A ~ above leads to 
some simplification Q. Thus, we have g oS ~ 2 and 



a. 



rr 



2a 



IR 



a 



RR) 



(15) 



where a- denotes the cross section of DiDj annihilation 
into possible SM final states. Since Dj R have the same 
interactions [Eq. (|T2"j) ]. we have a RR ~ a u for A ~ 0. 

The expression for a IT due to Higgs-exchange diagrams 
follows from its real darkon counterpart [l(| El j and so 



(o-j/Urei) = 



,2\2 



(16) 



where h is a virtual Higgs with the same couplings as the 
physical h, but with the invariant mass y/s = 2m D . and 
h — > Xi is any kinematically allowed decay mode of h. 
For m D > m h , the DiDj — > hh channel needs to be 
included [ll]. Thus the contributions of a T 
and J are 



'ii.rr t0 (<™) 



(°"/Z«rel) 



Jh 



(0"//«rel) 



2x, 



(17) 



For the Z'-mediated counterpart, assuming nonrelativis- 
tic darkons and A ~ 0, we derive 



a IR 



2 A™ 2 
D 



12n f (Am 2 D - m|,) 2 + r|,m|, 

where the sum is over fermions with mass < m D and 
N/ colors and Tz 1 is the Z' width. It follows that 



(18) 



36 



IR 



J z , 



36 



IR 



2 X ) 



(19) 



where b IR = a IR /v lel . Applying 

(av) = (av) h + (av) z , , J = J h + J z , (20) 

in Eq. (|14|) . we can then extract constraints on A^,^ and 
g ( from the 90%-CL range 0.092 < Vl D h 2 < 0.118 of the 
observed relic density [25j . 

In Fig. [21 we display as a function of m D subject to 
the relic data, assuming the absence of the Higgs effect 
and, for definitcness, m z , = 300 GcV, and compare it 
with the collider bound. This illustrates that only the 
resonance region, m z , ~ 2m D , can satisfy both sets 
of constraints. Since the situation is unchanged in the 
presence of the Higgs contribution, hereafter we limit the 
Z' contribution to this resonance case. If the former is 
nonnegligible, we can consider various combinations of 
Jh and J Z i with different implications for the prediction 
in relation to DM direct searches. 



0.500 - 

0.100- 
0.050 - 



0.010 
0.005 



0.001 



50 



100 



150 200 
m D (GeV) 



300 400 



FIG. 2: Curve of versus m D subject to the relic abun- 
dance data for m z , = 300 GeV and no Higgs-mediated con- 
tribution. The green solid (red dotted) section is allowed (ex- 
cluded) by collider data. 



B. Direct detection 

The direct detection of DM is through the recoil of 
nuclei after it hits a nuclcon N . In our model, with 
nearly degenerate Di n, the involved interactions are 



D m) N 



D I{R) N and D I(R) N 



Dr(i)N via h and 



Z' exchanges, respectively, in the t channel. The result- 
ing spin-independent cross-section a DN needs to accom- 
modate these possibilities, but also take into account the 
fact that the DM local density is independent of the num- 
ber of DM components. Accordingly, we find 



\ 2 2 2 2 

A DH 9NNh VDN v H 



'DN 



9q ^dn 



(21) 



in the nonrelativistic limit, where g NNh is the Higgs- 
nucleon effective coupling, ijl dn = m D m N / (m D + m N ), 
the first term is equal in form to the cross section in 
the real darkon case [l(| [13 j an d for the second term we 
have used (N\u^ a u+d'y a d\N) = 3Nj a N with vanishing 
contributions from the other quarks [26J . 

Here we look at a couple of representative examples. In 
the first one, the Z' contribution considerably dominates 
J in Eq. ([20]) . and we choose Jz> = 999 J;, for definitc- 
ness. We illustrate the prediction for a DN as a function 
of the darkon mass with m h = 125 GeV and m z , = 2m D 
in Fig.[3ja). For this plot, \dh follows from J = lO 3 ^ 
plus the relic constraint, g^ varies in the range allowed 
by the collider data, v H = 246 GeV, and we have em- 
ployed the range 0.0011 < g NNh < 0.0032 for the Higgs- 
nucleon coupling [l3|, H3|- The lightly shaded (lighter 
orange) portions of the prediction curve indicate that in 
the m D > 100 GeV region the Z' effect also dominates 
a DN and is enhanced by a few orders of magnitude rel- 
ative to the Higgs contribution. Evidently, compared to 
the most recent data from the leading direct searches 
for WIMP DM, the prediction can largely escape the 
strictest bounds to date. However, future direct searches 
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FIG. 3: Darkon-nucleon scattering cross-section correspond- 
ing to (a) J = 10 3 A (orange regions) and (b) J = 10 A 
(red regions) for m h — 125 GeV and m z , — 2m D . The 
darker (lighter) portions of the predictions come from the 
contributions of h alone (both h and Z'). The predictions are 
compared to 90%-CL upper limits from CDMS, XENON10, 
and XENON100 as well as the 90%-CL signal (pur- 

ple) region suggested by CoGeNT, a (gray) patch compati- 
ble with DAMA modulation signal at the 5a level, and two 
2cr-confidence (blue) areas representing CRESST-II data [29| . 
Also plotted is the XENON1T projected sensitivity [U . 



\/2 (2) times the corresponding coupling (cross section) 
in the simplest darkon model. We note that in all these 
instances darkon masses larger than ~100GeV are still 
viable and will be probed by future measurements. 



C. Invisible decay of the Higgs boson 



For 



> 2m,D R ,Dr, the decays h — > Dr jDri will 
occur and. for non- negligible Xdh, can substantially en- 
hance the invisible decay rate of the Higgs. Their com- 
bined branching ratio is 



B{h -> DiD! + D R D R ) 



T/i^DiDj + Th^yD R D R 



(22) 



where the Higgs width T h = lf u +r h ^ DlDl + T h ^D R D R 
includes the SM Higgs width r^ M . We depict in Fig.|l 
the branching ratios for the same Jh.z 1 and m h choices 
made in Fig. [3] and A ~ 0. Clearly, unless Jh is very 
small relative to Jz>, these invisible channels tend to 
dominate the Higgs width. This is similar to the real 
darkon case [To|, EH EH , in which the Higgs can be hid- 
den from sight for sufficiently low m h > 2m D . The 
model can thus readily account for the possibility that 
the LHC does not detect any Higgs. On the other hand, 
if an SM-likc Higgs is observed at ~125GeV, the model 
can also accommodate it and still provides light DM with 
the right relic abundance via the Z ! contribution. In con- 
trast, the simplest real-darkon model with m D < m h /2 
would be disfavored by such a discovery [23| . More gen- 
erally, depending on the parameters, our model can ex- 
plain a partially hidden Higgs, while leaving the ratios of 
branching fractions to SM particles essentially the same. 
For 2mu I ,D R > Tn h , the Higgs decay pattern would be 
SM-like. 

The recently discovered particle of mass near 125 GeV 
has properties consistent with those of an SM Higgs, 
based on the LHC information reported so far [l| . How- 
ever, the present data still allow the new boson to have 



such as XENON1T would probe this parameter space of 
the model more stringently. 

In Fig.[3jb) , we show an example where the Higgs con- 
tribution is less suppressed, Jh = Jz 1 /9. In this case, 
most of the m D < 100 GeV range is ruled out by the 
null results of some of the direct searches. Nevertheless, 
the prediction (red areas) is partly consistent with the 
possible WIMP hints reported by CoGeNT (purple ar- 
eas) and CRESST-II (blue areas), although they conflict 
with the other experiments. 

If the darkon annihilation is dominated by the Higgs 
contribution instead, Jh ^> Jz>, the allowed parameter 
space for m D < 100 GeV would be further reduced com- 
pared to that in the second example. With J ~ Jh, the 
relic density would comprise approximately equal parts 
from Di n and consequently X DH (cjjjv) would be about 
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FIG. 4: Total branching ratios of h —¥ Di,rDi,r versus m D 
corresponding to J = WJh and 10 3 Jh for m h — 125 GeV. 
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a branching ratio into invisible particles of up to a few 
tens of percent [3l|. All this is compatible with one of 
the possibilities discussed above, where the Z' contribu- 
tion greatly dominates the relic density. Obviously, fu- 
ture measurements at the LHC and DM direct searches 
together will test such a scenario within our model. 

V. SUMMARY AND DISCUSSION 

Anticipating the upcoming results of the Higgs hunt 
at the LHC and direct searches for DM in various under- 
ground experiments, we have considered a simple model 
possessing only a small number of nonstandard parti- 
cles, including scalar DM — the darkon — which is stabi- 
lized by a Z2 symmetry naturally arising from the spon- 
taneous breaking of a gauged U(1)b-l group. The asso- 
ciated gauge boson, Z' , yields new effects on the DM relic 
abundance in the resonant case, m z , ~ 2m D . We use 
the measurements of e + e~ — > £ + £~ and pp — > £ + £~X 
(£ = /.t, t) and the relic density data to constrain the new 
gauge coupling. Subsequently, we explore different com- 
binations of the Higgs- and Z'-mediated contributions 
to darkon annihilation subject to the relic density con- 
straints and evaluate the corresponding darkon-nucleon 
scattering cross-section, cr DN , compared to direct search 
results. If the Z' effect dominates the darkon relic den- 
sity, a DN tends to evade the present limits and may come 
mainly from the Z'-mediated contribution, depending on 
the darkon mass. We also discuss some implications of 



the invisible decay of the Higgs if the darkon is sufficiently 
light, m D < m h /2. In that case, we find that the Higgs 
invisible decay branching-fraction is still significant for 
most of the allowed m D values even if the Higgs contri- 
bution to darkon annihilation is very small compared to 
the Z' contribution. This result is important because the 
invisible branching fraction of the Higgs-like particle re- 
cently observed at the LHC has now been estimated from 
current data to reach up to a few tens of percent, which 
implies, within the context of our model, the necessity 
of a highly dominant Z' effect on the DM relic density. 
In general, the parameter space of the model has enough 
room to accommodate various potential outcomes of the 
ongoing Higgs and DM direct searches, which will there- 
fore probe the model further. The LHC can offer addi- 
tional tests via processes that can produce the Z' and/or 
a darkon pair, the latter due to diagrams mediated by 
h or Z 1 . The Z' may be most detectable in final states 
containing a pair of charged leptons, whereas the darkon 
may be uncovered in events with missing energy. 
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